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ABSTRACT 

In this paper we present the results of a high resolution (5") CARMA and SZA survey of the 3 mm 
continuum emission from 11 of the brightest (at 1.1 mm) starless cores in the Perseus molecular cloud. 
We detect 2 of the 11 cores, both of which are composed of single structures, and the median 3cr upper 
limit for the non-detections is 0.2 MqIu a '^5" beam. These results are consisent with, and as stringent 
as, the low de t ection rate of compact 3 mm continuum emission in dense cores in Perseus reported 
bv lOlmi et al.l ()2005[ ). From the non-detection of multiple components in any of the eleven cores we 
conclude that starless core mass functions derived from bolometer maps at resolutions from 10"-30" 
(e.g. with MAMBO, SCUBA or Bolocam) are unlikely to be significantly biased by the blending of 
lower mass cores with small separations. These observations provide additional evidence that the 
majority of starless cores in Perseus have inner density profiles shallower than r~^ . 
Subject headings: stars: formation; ISM: dust, extinction 



1. INTRODUCTION 

Surveys of starless cores in nearby molecular clouds 
have shown that the starless core mass function 
(CMF) looks very similar to the stellar initial 
mass function (IMF), but translat ed towards higher 



masses by a factor of a f ew (e .g. iMotte et al.l 
Nutter fc Ward-ThompsonI l2007t lEnoch et al 



1998 : 



20081: 



Simpson et al.ll2008D . If each starless core collapses into 



a single star, then the apparent shift between the IMF 
and the CMF implies that the formation process is in- 
efficient and most of the mass in each core never makes 
it into any star. Models have shown that 25%-70% of 
the mass in a protostellar core will be converted into a 
star with the remaining material being ejected by out- 
flows (Matzner & McKee 200^ . A second reason for the 
shift between the CMF and the IMF is likely to be frag- 
mentation, in which a single starless core will eventually 
form multiple stars. Since many stars are found in bi- 
nary or higher-order systems (e.g. iDuauennov fc MavoH 
Il991f ). to some extent the fragmentation of a starless core 
into multiple bound objects is to be expected. However, 
lower mass stars are less frequen tly found in gravitation- 
ally bound systems (|Ladall2006[) . 

Although the shape of the IMF can be reproduced 
from the CMF with a variety of assumptions about 
the multiplicity, fragmentation, and efficiency of the 
process (Swift & Williams 200g ?), recent work by 
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I Goodwin et al.l (|2008f) has shown that models in which 
every starless core eventually forms multiple stars repro- 
duces the stellar IMF more successfully than models in 
which low-mass cores are more likely to form a single 
star than a binary. That each core fragments into only 
two or three components, and not many more, is sug- 
gested by the relatively small number of close binaries 
and high velocity single stars that would result from the 
dyna mical ejection of pr otostars from higher order sys- 
tems (| Goodwin fc Kroup a 2005). Therefore, high angu- 
lar resolution observations of starless cores might be ex- 
pected to show that what appears to be a single object 
at low resolution is in fact composed of multiple frag- 
ments, depending on whether or not the fragmentation 
takes place before an embedded protostar forms and on 
what physical scale fragmentation occurs. 

The binary fraction of young stellar objects and pro- 
tostars has been studied in detail, and it has been found 
that these objects are commonly found in multiple sys- 
tems. For example, near-infrared observations of Class- 
I and flat-spectrum objects in nearby molecular clouds 
show that in the range 300-2000 AU the multiplicity 
is '^18%, with a medi an separation between compan- 
ions of ~900AU (Hais ch et all l2004f ). Sub-arcsecond 
Berkeley-IUinois-Maryland Array (BIMA) observations 
of the 3 mm continuum emission from 20 Class and 
Class I protostars found that all of the embedded objects 
are in small groupings, with separations ranging from 
<100AU to >10000AU (Looneyetal. 2000). Recent 
simulations of fragmentation in magnetically supported 
dense cores typically start with small cores ('^6000 AU or 
24" at the distance of Perseus) that fragment into multi- 
ple protostellar s ystems with separations on the order of 
10-1000 AU fe.g. [Bosil200l [20071 : IPrice fc Batdl2007h . 

Observations of starless cores have not often found 
evidence for fragmentation on the scale of 5-10". In 
a survey of dense cores in the Perseus molecular cloud 
with the Ow ens Valley Radio O bservatory (OVRO) in- 
terferometer, lOlmi et all ()2005f ) found that five out of 
six cores had no compact continuum emission, and in 
the sixth core a single source was detected (at ~5" res- 
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olution, or 1200 AU in Perseus). Similarly. iHarvev et alj 
([2OO1I) mapped the 1.3 mm continuum emission from the 
starless core L694-2 with the IRAM Plateau de Bure 
Interferometer (PdBI) and BIMA and found that the 
PdBI observations resolved out the core entirely while 
the shorter baselines of the BIMA data can be adequately 
modelled by a single source, implying that no fragmen- 
tation has taken place. Other interferometric observa- 
tions of starless cores, such as S68NW, L1544 and L694-2 
show only single peaks in their integrated spectral line 
emission (Wilhams & Myers 1999; WiUiams et al. 1999, 
120061 ). However, BIMA mapping of the N2II+ emission 
from the starless core L183 shows that this core is com- 
posed of three components with a total mass of a few 
tenths of a solar mass that are disti nct both spatially 
and kinematically (Ki rk et al.l 120091) . IKirk et all ()2009f ) 
interpret this as evidence that L183 is fragmenting as it 
collapses and predict that this core will eventually form 
multiple protostars. 

Here we present the results of a search for substruc- 
ture in a sample of 11 starless cores in the Perseus 
molecular cloud. We use Combined Array for Research 
in Millimeter-wave Astronomy (CARMA) and Sunyaev- 
Zel'dovich Array (SZA) interferometric observations to 
study the 3 mm continuum emission at ^-^5" resolution 
with sensivity up to scales of 30"- 80". The source selec- 
tion and observations are discussed in Section [2l followed 
by an analysis (Section [3|), a discussion (Section |4l), and 
summary (Section [5]). 

2. OBSERVATIONS 

To search for substructure within starless cores, we 
mapped the 3 mm continuum emission with ^5" resolu- 
tion towards a sa mple of 1 1 cores draw n from the 1.1 mm 
Bolocam survey of lEnoch et al.l (|2008[ ). Here we describe 
the sample of cores and their observations. 

2.1. Source Selection 

We drew our sa mple of starless core s from the 1.1mm 
Bolocam surve v oflEnoch et al. (2006, 2008), which has 
31" resolution. lEnoch et al.l (2008 ) uses Spitzer near- to 
mid-infrared data to determine the starless/protostellar 
status of each millimeter core. We chose starless cores in 
the Perseus molecular cloud on the basis of having high 
peak fluxes (>200 mJy beam~^). There are 15 star- 
less cores in the Bolocam survey of Perseus with peak 
fluxes >200 mJy beam^^, of which we were only allo- 
cated enough time to observe 1 1 , chosen arbitrarily. The 
selected cores have a median mass of 1.47 Mq and a me- 
dian density o f 3.2 x 10^ cm~ ^, as derived from the Bolo- 
cam data (En och et al.ll2008l ). The 1.1 mm-derived prop- 
erties of the starless cores in this sample are presented in 
Table [H 

Protostars in lEnoch et al.l ()2008D are identified by the 
shape of their near-infrared to far-infrared spectral en- 
ergy distributions (SEDs), minimum flux at 24 /im and 
presence of 70 /im point sources not classified as galaxy 
candidates. In Enoch et al. (200_8i), if a protostar is found 
within one FWHM of the center of the 1.1 mm Bolocam 
core, then it is classified as protostellar, otherwise the 
Bolocam core is classified as starless. Searches for proto- 
stars in the Perseus molecular cloud hav e also be en per- 
formed bv ,J0rgensen et al.. (|2007, ). ,Hatchell et al.l ([2007,}, 



and iSadavov et al.l ()2010D . |j0rgensen et al.l (|2007D iden- 
tify protostellar cores from their mid-infrared colors and 
magnitudes, distance to the nearest 850 /x m SCUBA core 
and co ncentration of the 850 /im emission. iHatchell et al.l 
(12007!) identify protostellar cores by the presence of ^^CO 
outflows and the presence of near-infrared (2MASS) or 
mid-i nfrared (Spitzer) sour ces within an 850 /im SCUBA 
core. ISadavov et al.l ()2010D find protostars based on mid- 
infrared and far- infrared (Spitzer) magnitudes and colors 
cross-referenced with a list of 850 /tm SCUBA cores. A 
detailed comparison between these methods of classify- 
ing a core as e ither starless or protostellar is given in 
ISadavov et al.l ()2010D : here we note that the classifica- 
tions according to the various criteria are broadly con- 
sistent. 

All o f the cores in our sample are classified as star- 
less by lEnoch e t al. (2008) and non e of the s e core s are 
identified as protostellar by J0rgensen et a l.l (J2007I) , but 
Pe rbo45 and PerboSO ar e identified as Class protostars 
bv IHatchell et al.l (|2007D ai id Perbo44 and Perbo l07 are 
identified as protostars by ISadavov et al.l ()2010l) . More 
recent deep 70 /im observations of one starless core in 
our sample, Perbo58, indicate that it may contain a very 
weak internal luminosity source (Enoch et al., in prep). 
Thus, Perbo58 may not be truly prestellar, although even 
if a central protostar has formed, it is likely very young 
and not far evolved from the prestellar phase. We keep 
Perbo44, Perbo45, Perbo50, Perbo58 and Perbol07 in 
the starless sample based on the original selection cri- 
teria and note that out of the four surveys to identify 
starless cores and/or protostars in the Perseus molecular 
cloud none of our targets are classified as protostellar in 
more than one paper. 

2.2. CARMA and SZA 

Continuum observations in the 3 mm window were ob- 
tained with CARMA, a 15 element interferometer con- 
sisting of nine 6.1 meter antennas and six 10.4 meter 
antennas. The CARMA correlator records signals in 
three separate bands, each with an upper and lower side- 
band. We configured all three bands for maximum band- 
width (468 MHz with 15 channels per band) to observe 
continuum emission, providing a total continuum band- 
width of 2.8 GHz. The observations were centered around 
102 GHz, and range from 100 GHz to 104 GHz. The half- 
power beam width of the 10.4 m antennas is 66" at the 
observed frequencies. Single pointings towards all 11 
starless cores were taken in the CARMA D-array config- 
uration, with projected baselines that range from 11m 
to 150 m. Follow-up observations of six of the eleven 
cores were made in the CARMA E-array cofiguration 
with projected baselines that range from 8 m to 66 m, also 
with single pointings. Two cores (Perbo45 and Perbo58) 
were then re-observed in D and E-array configurations 
with 7-point mosaics to improve the signal-to-noise in 
the maps and enlarge the area of uniform gain to 66". 
A summary of the CARMA configurations, the synthe- 
sized beam size, the largest angular scale to which the 
observations are sensitive, and the noise in the maps are 
presented in Table [2] 

The SZA is an 8 element interferometer equipped with 
3 mm band (80-115 GHz) receivers. The instantaneous 
field of view of the SZA is given by the primary beam 
of the 3.5 m antennas, approximately 3.3' (FWHM) at 
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95 GHz. Six of the SZA antennas are arranged in a close- 
packed configuration with baselines ranging from 4.5 to 
11.5 m and two outer telescopes yield baselines of up to 
65 m, resulting in a synthesized beam size of 38". Sixteen 
500 MHz wide analog bands provide 8 GHz bandwidth, 
though for the observations in this study only 7 GHz 
bandwidth (from 92-98 GHz) was availa ble. A more de- 
tailed description of the SZA is given in iMuchovei et ahl 



The observing sequence for the CARMA observations 
was to integrate on a phase calibrator (3C84, 0237+288 
or 0336+323) for 4-5 minutes and science targets for 
14-28 minutes. In each set of observations 3C84 was 
observed for passband calibration, and observations of 
a planet (Mars or Uranus) or 3C84 were used for ab- 
solute flux calibration. Although 3C84 has a variable 
brightness, it is monitored on a weekly basis at CARMA. 
Based on the repeatability of the quasar fluxes, the esti- 
mated random uncertainty in the measured source fluxes 
is cr ~ 5%. Radio pointing was done at the beginning of 
each track and pointing constants were updated at least 
every two hours thereafter, using either radio or optical 
pointing routines. Calibration and imagin g were done 
using the MIRIAD data reduction package (jSault et al.l 
[l99l . 

The observing sequence for the SZA observations was 
to integrate on the phase calibrator, 3C84, for 5 min- 
utes and the starless core Perbo58 for 10 minutes. The 
observations of 3C84 were also used for bandpass and 
absolute flux calibration. Because the SZA data are 
at a somewhat lower frequency than the CARMA data 
(95 GHz vs. 102 GHz), we use the flux of 3C84 derived 
at the higher frequency to set the absolute flux calibra- 
tion of the SZA data before combining the CARMA and 
SZA maps. Calibration and imagi ng were done usin g the 
MIRIAD data reduction package (|Sault et al.lll995f ). 

The 3 mm continuum maps of Perbo45 and Perbo58 
are presented in Figure [T] The other starless cores were 
not detected. 

2.3. Bolometer Maps 

In addition to the 31" Bolocam map at 1100 /im 
([Enoch et al.l [2006), the Perseus molecular cloud has 
also been mapped with SCUBA at 45 and 850 iim. 
at 9" and 14" resoluti on, respectively (iHatchell et al.l 
[2005t iKirk et al.l [20061 ). We downloaded the sub- 
millimeter maps from th e SCUBA Legacy Catalogues 
([Di Francesco et al1l2008[ ). which includes all of the data 
taken with SCUBA and uses the most current calibration 
and reduction methods. The absolute calibration of the 
bolometer maps is good to 15%, 20% and 50% at 1100, 
850 and 450 //m, respectively. The 450, 850 and 1100 //m 
bolometer maps of Perbo45 and Perbo58 are overlayed 
on the 3 mm interferometric maps in Figure [2j 

3. ANALYSIS 

We detected, at the >5cr level, continuum emission 
from 2 out of 11 starless cores in our survey. Here we 
present an analysis of the properties of the detected cores 
and the implications of this detection rate. 

3.1. Core Properties 

Both of the detected cores in this survey (Perbo45 and 
Perbo58) are seen to contain one main component (see 



Figure [T|). A few lower significance (3-4(t) peaks are 
present in the 3 mm map of Perbo58 along an extended 
ridge seen in the bolometer maps. Note that the map of 
Perbo58 includes data taken with the SZA, so it is sensi- 
tive to angular scales smaller than ^^80". Given that the 
FWHM sizes of Perbo58 as measured by Bolocam (26") 
is smaller than this scale, we believe that we have not re- 
solved out any of the 3 mm emission from Perbo58, but 
given that the map of Perbo45 only includes CARMA 
D+E array data and not SZA data we are likely resolv- 
ing out some of the 3 mm emission from Perbo45. 

3.1.1. Mass 

We calculate the mass of the detected cores from the 
total 3 mm flux and upper limits to the non-detections 
using the equation: 



M = 



^3mm 



B^(Td)k3„ 



(1) 



where d is the distance to Perseus, 6*3111111 is the total 
fiux from each detection or three times the noise for the 
non-detections, Bi, is the Planck function at dust tem- 
perature To, and Ksmm = 0.00169 cm^g"-^ is the dust 
op acity, extrapolated frorn the opacity at 1300 /xm given 
in lOssenkopf fc Henningj ([1994D Table 1, column (5) for 
dust grains with thin ice mantles, assuming a gas-to-dust 
ratio of 100 and an emissivity given hy (3 — 2. This value 
of the emissivity spectra l index is c onsistent with the 
measured P in TM C-IC ([Schnee et al. 2010 ) and L1498 
([Shirlev et al.ll2005D . For con i pariso n with the core prop- 
erties derived in lEnoch et al.l ([2008D . we assume that the 
distance to Perseus is 250 pc and that the dust temper- 
ature is 10 K for all of the cores in this sample, which 
is a typical gas tem perature of starless cores in Perseus 
([Schnee et al.ll2009l) . The derived masses and upper lim- 
its are presented in Table [31 

The mass of Perbo45 derived from our 3 mm data is 
about a factor of 4 lower than the mass derived from 
the 1.1 mm Bolocam data. The mass of Perbo58 derived 
from our 3 mm observations is a factor of 3 higher than 
the mass derived from the 1.1mm Bolocam data. Four 
possibilities for explaining these discrepancies are 1) un- 
certainties in the dust properties 2) variations in the tem- 
perature within the cores, 3) errors in the measurement 
of the 1.1mm and/or 3 mm fluxes, and 4) in the case of 
Perbo45, we are likely to have resolved out some of the 
3 mm continuum emission. Note that the deconvolved 
size of Perbo45 in the Bolocam surveys is roughly 53", 
so the amount of flux resolved out by our CARMA ob- 
servations is likely to be substantial. 

The masses of Perbo58 derived from the 1.1mm and 
3 mm data could be brought into agreement if the emis- 
sivity spectral index were given by /3 ~ 0.5 — 1 instead of 
(3 — 2, as we assumed. Although the two best measure- 
ments of the emissivity spectral index in starless cores 
suggest that > 2 is a much better fit than /3 < 1 
([Shirlev et al.l 120051: ISchnee et all |2010[ ). it is possible 
that grain growth leads to a shallower emissivity spec- 
tral index at the high densities likely to be found at the 
centers of prestellar cores. For instance, the emissivity 
spectral index in a sample of nearby pro tostars has been 
measured to be in the range 0.25 — 1.5 ([Arce fc SargentI 
[2006.) ■ and values in this range might be plausible for 
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starless cores give n the recent evidence fo r grain growth 
in starless cores (jSteinacker et al.ll2009[ ). Because the 
masses in Table |3l arc derived from the 1.3 mm opacity 
(jOssenkopf fc Heim ing 1994) extrapolated out to 3 mm, 
an incorrect assumption about the emissivity spectral in- 
dex would lead to masses off by a factor of a few. If 
the emissivity spectral index in starless cores is closer to 
(3 — 1 than /3 = 2, then the expected fluxes at 3 mm 
as extrapolated from the 1.1mm Bolocam observations 
would be a factor of '^3 larger than we have assumed. In 
this case the non-detection of 3 mm continuum emission 
towards nine of the eleven cores in our sample becomes 
an even more significant result (see Section [3. 2p . 

Alternatively, if Perbo58 were significantly warmer at 
the center than its average temperature, the higher res- 
olution 3 mm observations would be more biased to- 
wards high masses and densities than the low resolu- 
tion 1.1mm data. It is generally observed that starless 
cores are colder at their centers due to self-shielding (e.g. 
iCrapsiet ahl 120071: iSchnee et alll2007[) . Heating due to 
contraction can result in a temperature s a few degrees 
highe r in the center of a starless core (jKeto fc Casellil 
120091) ■ but this on its own cannot account for the large 
difference in masses derived from the 1 . 1 mm and 3 mm 
continuum data. If there were a weak internal luminosity 
source in Perbo58 (as discussed in Section [^?T|) . then one 
would expect some central heating that could bring the 
1.1mm and 3mm-derived masses into agreement. 

A third possible cause for the disagreement between 
the 1.1 mm and 3 mm derived masses is the uncertain ab- 
solute calibration of the emission maps. The two cores 
detected by CARMA/SZA are both in the dense and 
bright cluster NGC1333, which results in artifacts in 
the bolometer maps. However, detailed testing of the 
accuracy of the recovery of 1.1 mm fiux from Bolocam 
has s hown that the quo ted fluxes are accurate to within 
15% (jEnoch et al.ll2006D . The absolute calibration of the 
3 mm continuum maps is accurate to within 20%. There- 
fore, uncertainties in the absolute calibration of the emis- 
sion maps alone cannot account for the factor of a few 
difference between the 1.1 mm and 3 mm derived masses. 

Due to uncertainties in the dust emissivity and tem- 
perature profiles of Perbo45 and Perbo58, we will defer a 
more detailed analysis of their masses until more data are 
available and only note that the 3 mm fluxes of Perbo45 
and Perbo58 are not what we expected based on their 
1.1mm fluxes as mapped by Bolocam. 

3.1.2. Density 

We calculate the density of the detected starless cores 
from the derived mass and effective radius using the 
equation: 

n = ^ (2) 

where Toff is the effective radius, fj, — 2.33 is the mean 
molecular weight per particle and mn is the mass of 
hydrogen. The effective radius is the geometric mean 
of the deconvolved semi-major and semi-minor axes, de- 
rived from a Gaussian fit to the fiux distribution. The 
density for each detected core is presented in Table [3] 

The densities of Perbo45 and Perbo 58 derived from 
the 3 mm continuum data in this survey are larger than 
the values derived from the 1.1mm Bolocam data. The 



densities derived from the 3 mm data (see Table [3]) are 
roughly a few xlO^ — 10^ cm"'^, while the mean den- 
sity in a 10 AU diameter aperature derived from the 
1.1 mm data are on the order of a few x 10^ cm"-^. Since 
the 3 mm observations are higher resolution than the 
1.1mm observations, it is not surprising that the 3 mm- 
derived densities are larger. Though high, densities in 
the 10^ — 10^ cm~^ range have been reported in L1544 
( Crapsi et al.1120 07) an d in several starless cores in Orion 
( Nutter &: Ward-Thompsonll2007[ ). The densities in Ta- 
ble [3] would be lower if the assumed temperature were 
higher, as might be expected if the detected cores har- 
bor low-luminosity protostars. 

A simple model of the density profile of a starless core 
can be given by: 



n(r) = 



no 



1 



(3) 



where no, ro and a are constants. This formulation with 
a = 2.5 is a good appro ximation for the den sity profile of 
a Bonnor-Ebert sphere (iTafalla et al.ll200^ and provides 
a qualitatively good fit to the observed density profiles 
of starless cores, which are fairly flat at small radii and 
are steeper at larger radii. 

We investigate whether or not the density profiles of 
the cores in Perbo45 and Perbo58 are consistent with 
Equation [3] First we measure the fiux profile in annuli 
centered on the positions of Perbo45 and Perbo58 given 
in Table [3l We then calculate the fiux profile that would 
be observed for a spherical core with the density profile 
given by Equation [3] assuming that the cores are isother- 
mal and have constant dust emissivities. We normalize 
the observed and model fluxes such that the flux at the 
center of the core is equal to unity, so the no term in 
Equation [3] is not meaningful. We show in Figure [3] that 
under the stated assumptions of the dust properties and 
geometries that a density profile like that given in Equa- 
tion [3] does predict flux profiles that are similar to those 
observed at 3 mm for values of rp in the range of 5"- 
8"and a in the range 2.5-3.5. 

The model density profile that we use to fit the ob- 
served flux proflles of Perbo45 and Perbo58 does not pro- 
vide a unique solution to what the true density profiles 
of these core are, and other models that have a flat inner 
profile that steepens at larger radii can be consistent with 
our observations. If one uses different assumptions about 
the temperature profile and allows the emissivity to vary 
with radius, it might even be possible to fit the 3 mm 
flux profile with a power-law density profile. What we 
have shown here is that the density profiles of Perbo45 
and Perbo58 are not necessarily different from typical 
starless cores, and the reason that they are detected in 
our survey while other cores are not is most likely due to 
their relatively large masses for their size. 

3.2. Detection Rate 

Of the 11 starless cores that were observed with 
CARMA, 2 were detected (Perbo45 and Perbo58) and 9 
were not. Neither of the detected starless cores are seen 
to break into multiple components. The non-detection of 
the majority of the starless cores is not due to low sensi- 
tivity (the median upper limit for the non-detections is 
0.2 M0 - see Section [3.1.1l for details), but instead is likely 
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the result of a lack of compact sub-structures within the 
cores. If all of the mass in each of the non-detected cores 
had been contained within a few compact structures, ev- 
ery core would have been detected in the CARMA ob- 
servations (see Section 13.1.11 for an explanation of how 
the mass upper limits in Table [3] were derived) . For ex- 
ample, the median Scr upper limit on the mass of the 
non-detected cores is a factor of 7 smaller than the core 
mass derived from the Bolocam map. Therefore, we ar- 
gue that the non-detection of multiple components in all 
11 cores in our sample is a significant result. Further- 
more, given the hi gh rate of multiplic ity in Class and 
Class I protostars (jLoonev et al.|[200(). ) and the low rate 
of multiplicity in starless cores, we suggest that fragmen- 
tation into binaries takes place during the collapse of a 
prestellar core or during the Class stage. 

The density profiles of many starless cores have 
been measured from (sub)millimeter continuum emission 
maps and it has generally been found that the density 
profiles are consistent with an r~^ power-law at large 
radii that flattens towards the center of the core (e.; 
IWard-ThomDson et all 11991 [T999t IShirlev et al.ll200f 
To study the density profiles of the starless cores in our 
sample, we use the mass and radius derived from the 
1.1 mm observations to predict the 3 mm fiux at the cen- 
ter of each core with the assumptions that the density 
profile goes as r~^, that the emissivity spectral index is 
given by /3 = 2, and that the cores are isothermal. We 
would have been able to detect (at the 3cr level) all 11 
cores in our sample if these assumptions were valid, but 
in fact we only detect two starless cores. As noted in 
Section I3.1.1[ if /3 = 1 then the detection rate of cores, 
and any substructure in them, would be easier than if 
/3 = 2. Because we expect starless cores to be colder in 
their centers than at the edges the flux profile of a "real" 
starless core will be somewhat flatter than our isothermal 
model would predict. Nevertheless, we conclude that the 
density profiles for the majority of the starless cores in 
our sample are flatter than an r~^ proflle at their centers, 
in agreement with other surveys. 

The low rate of 3 mm continuum detections towards 
11 of the brightest (>200 m Jy beam~^ at 1.1mm with 
Bolocam: lEnoch et al.ll2008D in this st udy is in agreement 
with the results of lOlmi et al.l ()2005[ ) who detect only 1 
out of 6 dense cores in Perseus. The largest angular 
scale to which the interferometer was sensitive and the 
median noise in the maps in this study (30"-40" and 0.7 
mJy beam^^) are similar to that in lOlmi et al.l (20"-40" 
and 1-1.5 mJy beam~^: l2005() . There is no ov erlap in the 
sampl e of sources mapped in this study and in lOlmi et al.l 
(|2005f ). 

There are 67 starless cores detected by lEnoch et al.l 
(J2008) in the Perseus molecular cloud, of which we ob- 
served 11. Because we selected our sample based on their 
high peak fluxes, there is no reason to expect that the 
other starless cores in Perseus would have a higher detec- 
tion rate at 3 mm. The low rate (0/11) of finding mutiple 
3 mm components at 5" resolution within 1.1mm cores 
observed with 30" resolution is consistent with the 3 mm 
continuum OV RO observa. t ions o f dense cores in Perseus 
carried out by lOlmi et all ()2005D , who found that none 
of the cores (out of six) were seen to be composed of 
multiple components. We argue that starless core mass 



functions ge nerated from single dish bolometer maps o f 
Perseus (e.g. lEnoch eFaLll2008l: iHatcheh fc Fulleiil2008D 
are not strongly biased by the blending of compact, close 
cores into more massive single cores. 

4. DISCUSSION 

Giv en that many stars are f ound in multiple sys- 
tems (|Duauennov k, Mavoiill99il) and that many proto- 
stars are also multiples (jLoonev et al.ll2000l : iHaisch et al.l 
|2004| ). it is worthwhile to examine when fragmentation 
takes place. If starless cores fragment before formation of 
an embedded object, then we would expect to be able to 
detect these sub-cores. On the other hand, if fragmen- 
tation takes place after the first protostar has formed, 
then a survey of starless cores won't find any evidence 
of fragmentation. In this study, we find that none of 
the eleven starless cores show any evidence of fragmen- 
tation. If these cores will eventually fragment, then this 
must happen during and/or after the collapse and forma- 
tion of the first protostar. Late fragmentation on small 
scales, which is consistent with the results of this survey, 
is predicted by theories that create binaries fro m disk 
fragm entation (e.g. Bonnell & Bate 1994; Kratte r et al.l 
l2010h . 

The observations described in this paper are in agree- 
ment with the re sults of other recent surveys in Perseus. 
iKirk et al.l ()2007[ ) compared the line profiles of C^^^O (2- 
1) and N2H"'"(l-0) toward a sample of 150 candidate 
dense cores, both prestellar and protostellar. They found 
that the N2II+ lines, representing the densest gas within 
the core, display nearly thermal line widths, indicating 
little non-thermal activity at the core center. While the 
C^®0 gas shows significant non-thermal broadening, the 
velocity centroids of the two species are very similar, in- 
dicating that the connection between the cores and their 
environment is relatively calm. Interestingly, such qui- 
escence across scales appears hard to r econcile with nu- 
merical simulations of turbulent clouds (|Kirk et al.ll2009D 
and may indicate a need for magnetic pressur e sup- 
port w ithin the dense gas. Additionally, Jorgenscn et al.l 
(|2007f l compared the location of mid-infrared Spitzer 
sources and sub-millimeter cores in Perseus, finding ex- 
cellent agreement between the deeply embedded proto- 
stars and their proto-stellar envelopes. Intriguingly, only 
3 out of 40 cores were found to harbour multiple proto- 
stars, at the spatial resolution of Spitzer (MIPS 24 /im, 
6" or '^ISOOAU). Thus, it appears that within Perseus 
there is little evidence for fragmentation within cores, 
at least until the protostellar stage, where the fragmen- 
tation length scale is smaller. This may be due to the 
smooth and quiescent manner in which the cores them- 
selves appear to be forming. 

We only detected significant (>5cr) emission from two 
of the eleven starless cores in our sample and we put 
tight upper limits on the masses of multiple components 
in the other nine cores, from which we conclude that 
there is no substructure on ~5" scales (~1200AU) in 
the majority of starless cores. A smooth column density 
profile is predicted by many models of isolated star for- 
mation, such as the insi de-out collapse model ( Shu 19771) 
and ambipolar diffusion ([Ciolek fc Mouschovia£Ll993.) . It 
is less clear if a smooth column density profile is con- 
sistent with cores created through turbulence, though 
some simulations show that turbulent molecular clouds 
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do form quiescent cores with smooth column density 
profil es (jNakamura &: Li|[2005l : iBallesteros-Paredes et alj 
[2OOI . 

5. SUMMARY 

We present a study of the 3 mm continuum emission 
from a sample of 11 of the brightest (at 1.1 mm) starless 
cores in the Perseus molecular cloud. The data were 
primarily taken with CARMA at ~5" resolution, with 
followup observations of one core taken with the SZA to 
improve the sensitivity to extended flux. 

1) We detect two starless cores (Perbo45 and Perbo58), 
both of which are seen to contain a single, main compo- 
nent. 

2) Based on the low detectio n rate (0/11 in this study 
and 0/6 in lOlmi et"aLl I2005D of multiple components 
within starless cores, we conclude that starless core mass 
functions derived fro m bolometer maps of the Perseus 
mole cular cloud (e.g. lEnoch et al.l l2008t iHatchell et al.l 
l2007f) are unlikely to be heavily biased by the blending 
of smaller, nearby cores into single sources at the coarser 
resolutions of single-dish maps. 

3) Based on the low detection rate of any compo- 
nents within den se cores (2/11 in this study and 1/6 in 
lOlmi et ani2005[ ). we conclude that the majority of star- 
less cores lack compact structures. More specifically, we 
find that density profiles in the inner portions of starless 
cores of the form r~^ are inconsistent with our oberva- 
tions for the majority of the starless cores in our survey. 
This agrees with the results of previous surveys, which 
find that the density profiles of starless cores are flatter 



than r ^ at small radii (e .g. lWard-Thompson et al.lll994 
[T999I: IShiriev et al.l[200(l . 
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Substructure In Starless Cores 
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Fig. 1. — 3mm flux maps (grayscale and contours) of Perbo45 (top)and Perbo58 (bottom), from data taken with CARMA and, for 
Perbo58, SZA. Contours start at 3(7 and increase by la, negative peaks are shown with dashed contours. The beam size is shown in the 
bottom left corner of each map. The units of the maps are mjy beam~^. 



TABLE 1 
Core Properties from Bolocam 1 mm Data 



Name 



RA 
(J2000) 



Dec 
(J2000) 



Peak Flux 
(mjy/beam) 



Total Flux 

(Jy) 



Total Mass 
(Mq) 



Derived from an elliptical Gaussian fit, deconvolved by the 31 beam 
Mean density calculated in a fixed linear aperture of diameter 10 AU 



FWHMi 

n 



Density^ 
(10^ cm-3) 



Perboll 


03:25:46.0 


+30:44:10 


241 


0.90 


2.18 


72.5 


2.2 


Perbol3 


03:25:48.8 


+30:42:24 


407 


0.47 


1.14 


38.2 


3.9 


Perbol4 


03:25:50.6 


+30:42:01 


342 


0.41 


1.00 


38.7 


3.5 


Perbo44 


03:29:04.5 


+31:18:42 


274 


0.72 


1.73 


47.0 


2.7 


Perbo45 


03:29:07.7 


+31:17:17 


455 


1.35 


3.25 


52.9 


4.8 


Perbo50 


03:29:14.5 


+31:20:30 


313 


1.31 


3.15 


68.0 


3.2 


PerboSl 


03:29:17.0 


+31:12:26 


423 


0.62 


1.49 


34.7 


3.6 


Perbo58 


03:29:25.7 


+31:28:16 


273 


0.33 


0.78 


25.9 


1.8 


Perbo74 


03:33:01.9 


+31:04:32 


255 


0.38 


0.91 


44.5 


2.6 


PerbolOS 


03:43:57.8 


+32:04:06 


283 


0.61 


1.47 


46.6 


2.8 


Perbol07 


03:44:02.1 


+32:02:34 


388 


0.48 


1.17 


46.0 


4.1 
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Fig. 2. — 3inm flux maps (grayscale) of Perbo45 (top) and Perbo58 (bottom), with the same scale as in Figure[T] Contours show the 
emission from maps made by Bolocam 1.1mm (left), SCUBA 850 /^m (middle) and SCUBA 450 /xm (right), with contour levels at 30%, 
50%, 70%, and 90% of the peak in each map. 
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Fig. 3. — 3 mm flux profiles of Perbo45 (left) and Perbo58 (right), shown in solid black lines. The flux profiles of model cores with density 
profiles given by Equation |3] are shown in dotted grey lines. The inner point on the observed mass profiles is placed at 1.25", which is half 
of the HWHM of the typical beam size of the 3 mm observations. Values of rg and a, as described in Section l3.1.2l are given in each panel. 
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TABLE 2 
3 mm Observations 



Name 


CARMA D 


CARMA E 


SZA 


Beam Size^ 
(") 


LAS2 
(") 


Noise 
(mjy/beam) 


Perboll 


Y 


N 


N 


4.5 


30 


0.5 


Perbol3 


Y 


Y 


N 


6.6 


40 


1.3 


Perbol4 


Y 


Y 


N 


6.2 


40 


0.9 


Perbo44 


Y 


N 


N 


4.2 


30 


0.6 


Perbo45 


Y 


Y 


N 


6.4 


40 


0.4 


Perbo50 


Y 


N 


N 


4.2 


30 


0.7 


Perbo5l3 


Y 


Y 


N 


5.4 


40 


2.8 


Perbo58 


Y 


Y 


Y 


5.6 


80 


0.3 


Perbo74 


Y 


N 


N 


4.6 


30 


0.3 


Perbol05 


Y 


N 


N 


4.3 


30 


0.9 


Perbol07 


Y 


Y 


N 


6.2 


40 


1.1 



Geometric mean of the major and minor axes 
Largest angular scale to which the map is sensitive 

The noise in the map of Perbo51 is larger than the other maps due to the bright source 
NGC 1333 IRAS4 located at the edge of the CARMA map 



TABLE 3 
Gaussian Fits to 3 mm Observations 



Name 


RA oflFseti 


Dec offset^ 


Peak Flux2 


Total Flux^ 


Axes'^ 


OPA^ 


Mass^ 


density 




(") 


{") 


(mJy/beam) 


(mjy) 


(") 


(degrees) 


(Mo) 


(em-3) 


Perboll 














<0.11 




Perbol3 














<0.29 




Perbol4 














<0.20 




Perbo44 














<0.14 




Perbo455 


-10.3 ±0.5 


8.8 ±0.7 


2.4 ±0.3 


11 ±0.5 


14 X 9 


-14 


0.8 


1.1 X lO'^ 


Perbo50 














<0.16 




Perbo51 














<0.62 




Perbo58 


-4.3 ±0.8 


-1.1 ±0.9 


2.0 ±0.3 


33 ±1 


26 X 18 


35 


2.4 


4.5 X 10" 


Perbo74 














<0.07 




Perbol05 














<0.20 




Perbol07 














<0.24 





Offset measured from pointing center, 



given m 



Table [T] 



Errors derived from the random noise in the maps and do not include absolute flux calibration uncertainties 
Deconvolved major and minor axes and position angle 
For non-detections, the 3(T upper limit to a point source mass is given 

Note that only CARMA D and E-array obserations have been taken for this source, so we are likely to have resolved 
out a larger fraction of the emission than for Perbo58 



